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Abstract 

We present a simple method to deal with caustics in the semiclassical 
approximation to the partition function of a one-dimensional quantum 
system. The procedure, which makes use of complex trajectories, is ap- 
plied to the quartic double-well potential. 



1 Introduction 

It is well known |Q that the partition function for a particle of mass m interacting 
with a potential V(x) and a thermal reservoir at temperature T can be written 
as a path integral ((3 = l/fceT): 

Z(p)= [ dx {x Q \e-^\x ), (1) 



(x \e-P H \x ) = / [Dx(r)]e- S ^ h , (2) 

Jx(0)=x 



S[x] = / cIt 



o 



— mx 2 + V(x) 



(3) 
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The semiclassical approximation to the density matrix element is given byg 

N 

(zo|e-^|zo>«£e- s ^ A -V 2; (4) 
fe=i 

where x^(t) is a classical trajectory [i.e., it is a solution to the Euler-Lagrange 
equation, mx = V'(x), subject to the boundary conditions x(0) = x((3h) = xq] 
that minimized^ (globally or locally) the action S[x], and is the determinant 
of of the fluctuation operator F[x^] = —md^. + V"[x*]. (A derivation of this 
result will be sketched in Section ^.) 

In Rcf. gj we have examined, for the sake of simplicity, potentials for which 
there is only one classical trajectory satisfying the above boundary conditions. 
In general, however, the number N of such solutions depends on xq and (3. A 
problem then occurs when we cross a caustic [the frontier between two regions 
of the (xo, j3) -plane characterized by different values of N]\ the r.h.s. of (|]) 
diverges jjjj. This divergence, however, is unphysical, being an artifact of the 
semiclassical approximation. The purpose of this work is to present a simple 
extension of the semiclassical approximation which circumvents this problem.^ 
(Due to limitations of space, here we shall only sketch the method. Details will 
be given elsewhere.) 



2 Improved semiclassical approximation 

In order to show how one can improve the semiclassical approximation so as to 
eliminate the unphysical divergences at the caustics, it is convenient to recall 
how (|J) is derived. Briefly, one has to: (i) expand the action around a minimum 
x c (t): S[x c + 77] = S[x c ] + S 2 + SS, where S 2 = \ f£ H dr V (t) F[x c (t)} V (t) and 
5S = 0(rj 3 ); (ii) throw away 8S; (iii) express t){t) in terms of the orthonormal 
modes of F, i.e., t](t) = Y^jLo a j t fj{ T )i where Fipj(r) = Xj<fj{r), ^j(O) = 

Vj(pK) = 0; then S 2 = \ E 3 °l X i a ) and [ Vx ( T )\ = WjLo da o/^^- The P ath 
integral in (||) has now become a product of Gaussian integrals. Performing 
the integrations one arrives at the "usual" semiclassical approximation to the 
density matrix element: 

(x \e-e»\x )*e- s ^A-^, (5) 

where A = YlJLo = det F. If there are N minima, one has to add together 
their contributions, thus obtaining (^). 

When we cross a caustic, a classical trajectory x c (t) is created or annihilated. 
Precisely at this point, the lowest eigenvalue of F[x c ] vanishes, thus making the 
integral J_ dag exp(— Aoa§/2?i) diverge. This problem can be remedied by 

x Each term in the sum on the r.h.s. of (^) is in fact the first term of a series. See Ref. [||. 
2 The Euclidean nature of the path integral allows one to discard saddle-points. 
3 Anker hold et at M] have discussed the caustics problem near the top of a potential barrier. 
The present work shows how to deal with the problem everywhere. 
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retaining fluctuations beyond quadratic in the subspace spanned by ipo (the 
eigenmode of F associated with Ao)- As a result of this procedure,^ we obtain 
an improved approximation to the density matrix element (|^): 

(x \e-^\x ) « e - s ^ h A- 1 / 2 F(x , /?), (6) 
where x gm (r) is the global minimum of S[x] and 



with 



F( Xo ,P) = x /^L I daoe-v^)/*, (7) 



V(oo) = \ Aoa 2 , + > J ( / r/r\-""[,, m (r)] -(>.») ^f. (M 




A couple of remarks are in order here: (i) we take for M [Eq. (jg)] the smallest 
even integer such that the coefficient of a^f in V(ao) is positive for all values of 
Xq and /?; this suffices to make the integral in (pp finite even when Ao vanishes; 
(ii) the factor Ag^ 2 in T cancels the factor A contained in A -1 / 2 ; combined 
with (i), this shows that the improved approximation to (xo\e~^ H |xo), Eq. (|J), 
is finite at the caustics; (iii) there is a one-to-one correspondence between the 
minima of S[x] and the minima of V(ao); therefore, it is not necessary to ex- 
plicitly add their contributions as in (Q), for they are already contained in T . 

Although the procedure outlined above teaches us how to cross the caustics, 
it is not very convenient: in order to obtain the coefficients of V(ao) one has to 
find Ao and ipa(r). This, in general, is not an easy task, and makes the whole 
procedure very cumbersome. Instead, we shall present an alternative way of 
obtaining those coefficients, which is based on remark (iii) above. 

Let us assume that M — 4 in Eq. (||); this is the case for the quartic double- 
well potential, to be discussed in the next section. Then the "effective action" 
A(ao) = S[xgm] + V(do) for the "critical" mode cpo is a fourth degree polynomial 
in ao. Let us assume it has three (real) extrema: a global minimum at ao = 0, 
a local maximum at u > 0, and a local minimum at v > u. This allows us to 
write such a polynomial as^| A(ao) = S[x gm ] + a [5 uv a 2 , — | (u + v) aj) + \ clq] ■ 
In order to relate the parameters in A(ciq) to calculable quantities, we impose 
that A(v) = 5[xi m ] and A(u) — S[x sp ], where xi m (r) and x sp (r) are the local 
minimum and the lowest saddle-point of S[x], respectively. This yields 

5[x lm ] - S[x em ] = A(v) - A(0) = g 3 (2-g) 

S[x sp ] - S[x gm ] A{u) - A(0) 2^-1 ' U 

where £ = v/u. Another combination of parameters which becomes determined 
is /i = an 4 : 

S[x sp ] - 5[x gm ] = A(u) - A(0) = ii (2g - 1). (10) 

4 The procedure adopted here resembles the treatment of caustics in optics |g| and in 
quantum mechanics ^j, suitably modified to take into account the Euclidean nature of the 
path integral in (Eh. 

5 One can easily check that A'(0) = A'(u) = A'(v) = 0. 
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In terms of the calculable parameters £ and fi, we can rewrite A(ao) as A(ao) = 
S[x gm ] + V 3 (ao/u), where V 3 (z) = p [§ £z 2 - § (1 + £) z 3 + \ z 4 ] . There still 
remains one parameter to be determined, namely u. Fortunately, T does not 
depend on it. Indeed, identifying V3(ao/u) with V(ao) yields A = /i£/u 2 , which 
allows us to rewrite (fij) solely in terms of £ and fx: 

V 2nTiJ_ 00 

The discussion above can be easily adapted to the case in which ^4(ao) has 
only one extremum. In this case, A'(ao) has one real (ao = 0) and two complex 
conjugate roots (w and w*), the latter corresponding to the complex trajec- 
tories x c t(r) and x* t (r). Accordingly, one has .A(ao) = S[x gm ] + Vi(ao/\w\), 
where Vi(z) = x [§ 2:2 — § (cos^) z 3 + ^ z 4 ] , with x = a\w\ 4 and <f) = arg(u>). 
Identifying A(w) with 5[x c t] then yields 

S[x ct ] - S[x gm ] = A(w) - A(0) = ^ (2w* = £ (2e 2 '^ - e 4 '^) , (12) 

from which we can obtain \ an d <j>- Finally, identifying Vi(ao/|u>|) with V(ao) 
leads to A = x/|w| 2 , so that 

f — I X [°° rf 0n# .-Vi(ao/h|)/B _ /_X_ [°° fl zp -Vx(z)/% do) 



' W 2^?m 2 



3 Application: the quartic double-well potential 

Let us consider the quartic double-well potential, V(z) = ^ (a; 2 — a 2 ) 2 , A > 0. 
In order to simplify notation, we replace x and r by q = x/a and = lot, 
respectively, where ui = (Aa 2 /™) 1 / 2 . In the new variables, the equation of 
motion reads q = U'(q), where U(q) — \ (q 2 — l) 2 . Closed trajectories have the 
form 

q c (9) = q t cd(u,k), (14) 

where cd is one of the Jacobian elliptic functions [Q, u = \]\ — q\ j2(Q — 0/2), 
= /3hu), and k = yjq 2 /{2 — q 2 ). The turning point qt is hxed by the boundary 
condition g c (0) = go- The classical action can be written as S[x c ] — (h/g) I[q c ], 
where g = Ti\/m 2 L0 3 and 

I[q c ] = 9 U(q t ) + 2sgn(g t - q ) /"* ^2 [17(g) - dg. (15) 

Finally, the determinant of the fluctuation operator is given by [p| 



A = 4 W (g - g t ) _ & . (16) 

Using these ingredients one can compute both the usual and the improved 
scmiclassical approximations to (go|e _/3H |go), Eqs. (0) and (||), respectively. The 
results are compared in Fig. 1. 
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Figure 1: p = (qo\e~^ H \qo) vs. qo for 6 = 5.0 and g = 0.3. The dashed 
line is obtained using the usual semiclassical approximation [Eq. (Q)]; the solid 
line is the result of the approximation discussed in Section ||. The equation 
g c (0) = qo (which determines the turning point qt) has three real solutions for 
Qo < <7o — 0.33319 . . .; as we cross the caustic, two of them (the ones associated 
with a local minimum and a saddle-point of the action) coalesce — thus causing 
the divergence in the usual semiclassical approximation — and reemerge at the 
other side as a pair of complex conjugate solutions. 
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